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We investigate the Mott transition using a cluster extension of dynamical mean field theory 
(DMFT). In the absence of frustration we find no evidence for a finite temperature Mott transition. 
Instead, in a frustrated model, we observe signatures of a finite temperature Mott critical point 
in agreement with experimental studies of k— organics and with single site DMFT. As the Mott 
transition is approached, a clear momentum dependence of the electron lifetime develops on the 
Fermi surface with the formation of cold regions along the diagonal direction of the Brillouin zone. 
Furthermore the variation of the effective mass is no longer equal to the inverse of the quasi particle 
residue, as in DMFT, and is reduced approaching the Mott transition. 



The Mott transition is one of the central issues in the 
physics of strongly correlated electron systems. Great 
theoretical progress in this area has been achieved using 
DMFT 0. In spite of these successes, the single site 
DMFT approach has several limitations due to the k- 
independence of the self-energy: variations of the quasi- 
particle residue, the quasiparticle lifetime and the effec- 
tive mass along the Fermi surface are absent. Moreover, 
some orders involving several sites like d-wave supercon- 
ductivity can not be described. Furthermore, single site 
DMFT does not capture the effects of the magnetic ex- 
change interaction on the single particle properties in 
the paramagnetic phase. This is potentially a singular 
perturbation, which can bring substantially new effects. 
Hence the question of whether the Mott transition sur- 
vives beyond the single site mean field approximation and 
how it is modified by short range magnetic correlations 
is among the most pressing questions in the field. 
These questions can now be addressed using cluster ex- 
tensions of DMFT III this paper, we report a study 
of the Cellular DMFT (CDMFT) approximation on 
a 2 X 2 square cluster of the one-band Hubbard model 
on a square lattice with and without frustration, solving 
the CDMFT equations with the Quantum Monte Carlo 
(QMC) method. In the unfrustrated model, the finite 
temperature Mott transition is precluded by antiferro- 
magnetic (AF) order. However, if the frustration is big 
enough to destroy this order, we show the existence of a 
Mott transition which is in agreement with the DMFT 
scenario and displays, at the same time, new physical 
effects due to the fc-dependence of the self-energy in the 
cluster method. In particular, we observe the appearance 
of different behavior of the one electron spectral function 
in different regions of the Fermi surface. This behav- 
ior, out of reach of single site DMFT, is very important 
since observed in the photoemission experiments 0, 01 of 
cuprate superconductors. 

The one-band Hubbard model on a two-dimensional 



square lattice with hopping tij 
U reads : 



and Hubbard repulsion 



«a-2) (1) 



and 



denote 



where a =1, i is the spin index, and ^ 
the electron operators. The bare one electron dispersion 
is given by Ej^ = —2t{cos{kx)+cos{ky)) — 2t'coa{kx + ky). 
We analyze both the unfrustrated model i' = and the 
strongly frustrated one f « t, where the antiferromag- 
netic J_AF) order is expected to be destroyed by frustra- 
tion We take D = 4t ~ 1. This type of anisotropic 
band structure is particularly relevant for the study of 
organic materials hke k-(BEDT-TTF)2X (see e.g. 0). 
Within CDMFT we divide the lattice in a superlattice 
of two by two squares and we basically solve the DMFT 
equations applied to the superlattice (See Fig. 
Starting from an effective action containing a Weiss dy- 
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FIG. I: The real space cluster and labeling conventions. 

namical field Gq one computes by QMC (Hirsch-Fye algo- 
rithm) a local cluster propagator and a cluster self-energy 
(which are four by four matrices). The Weiss field is then 
recomputed using the Dyson equation 

\ -1 
1 \ 



+ /i — i{K) — T,{iujr} 



S(ia;„) 



where K is in the reduced Brillouin zone of the super- 
lattice, t is the hopping matrix for the superlattice and 
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jj, is the chemical potential. This procedure is iterated 
until convergence is reached. In our QMC runs we use 
30000 — 300000 QMC sweeps and a number of time slices 
L « PU. Finally, let us recall that in CDMFT, the lattice 
self-energy is different from the cluster self-energy E and 
reads here {k is in the Brillouin Zone of the square 
lattice) : 



4^ " ' 2 

i=l 



(E12 + 1134) cos{kx 



(i;24 + Si3) C0s(fcy)-|-Ei4 COs{kx + ky)+Y.23 C08{kx-ky) 



Let us first focus on the unfrustrated case {t' = 0) at 
half-fiUing (See also [ill[lS[Ili). We computed the Neel 
temperature, corresponding to the continuous onset of an 
antiferromagnetic solution, as a function of U. Our result 
is similar, even quantitatively, to the one obtained in a 
DMFT study of the bipartite half-filled Bethe lattice [H. 
In particular the AF transition takes place at a quite high 
temperature, precluding a finite temperature DMFT-like 
Mott transition {i.e. a first order transition ending at a 
finite temperature critical point). 

Let us emphasize a major difference between single- 
site and cluster methods. In DMFT, the paramagnetic 
equations {i.e. the equations with enforced spin symme- 
try) are the same for the bipartite half-filled Bethe lattice 
and for a fully frustrated infinite dimensional model Q. 
Therefore the study of the Mott transition can be per- 
formed analyzing the paramagnetic solution in the AF 
part of the phase diagram. On the contrary, in cluster 
methods, the paramagnetic equations are not the equa- 
tions of such a simple fully frustrated model. Because of 
the quantum short-range fluctuations existing inside the 
cluster, studying the paramagnetic solution for parame- 
ters where an ordered solution exists is unphysical. In- 
deed if one forces the symmetry and follow the paramag- 
netic CDMFT solution inside the AF region a pseudo-gap 
opens because of the Slater mechanism. Some previous 
work have focused on increasing the cluster size H, . 
The Neel temperature then decreases to 0, but the T = 
AF order can still induce long range fluctuations at finite 
temperature and open a pseudo-gap. 

On the contrary, in this paper we analyze a strongly 
frustrated Hubbard model with a 2 x 2 cluster. In the 
following, we will take t'/t = 1. This frustration is 
big enough to destroy the AF order, at least for inter- 
mediate temperatures T/D > 0.02 where we have per- 
formed QMC computations. In fact we observe unam- 
biguous signatures of a finite temperature Mott transi- 
tion similar to the single site DMFT case. In Fig. |21we 
present the imaginary parts of the on-site cluster self- 
energy and propagator H'd and G'd, for different values 
of U and as a function of the Matsubara frequencies. For 
U < 2.2 and small w„, we observe a metallic behavior: 




FIG. 2: Imaginary part of the on-site propagator and of the 
on-site self-energy (inset) in Matsubara frequencies for U = 
2,2.1,2.2 (green curves with triangles), U = 2.25,2.29,2.31 
(blue curves with squares), U = 2.34,2.43,2.5 (red curves 
with circles) and T/D = 1/44. 



E"j(iw„) ~ (1 — -i) a;„. When w„ ^ 0, the imaginary 
part of the on-site self-energy increases linearly, while 
G'li decreases to a constant. This is the behavior charac- 
teristic of a Fermi liquid with a decreasing quasi-particle 
residue Z and a finite density of states at zero frequency. 
For U > 2.34 the behavior is clearly different: when ojn 
decreases to 0, E'd decreases while G'd increases. This 
is characteristic of an insulating phase with a (pseudo-) 
gap in the density of states. These results shows clearly 
the crossover from a metallic to an insulating regime, 
similarly to DMFT JJ. Furthermore, we have found in- 
dications of the existence of a finite temperature Mott 
critical point, located around Uc/D ~ 2.3 — 2.4 and 
Tc/D ~ 1/44. In Fig. El we present the double occupa- 
tion docc = \ Tlii=o ^ function of U for various 
temperatures. It presents a singular behavior similar to 
the one found in [l2| for DMFT. Furthermore the solution 
of the CDMFT equations displays a numerical "critical 
slowing down" around Uc,Tc. This is natural around a 
second order transition where new solutions appears con- 
tinuously. 

Let us now turn to the new aspects of the Mott tran- 
sition that one can get with cluster methods and that 
are excluded from the beginning within a DMFT analy- 
sis. As discussed previously, they are encoded into the 
/c— dependence of the lattice self-energy. First, let us fo- 
cus on the Fermi surface that can be located analyzing 
the maxima, as a function of fc, of the spectral weight 
at zero frequency A{k,u} = 0) « ^ lim^^^o ImG'(fc, zu;„). 
Since the QMC method produces finite temperature data 
in imaginary time, the self-energy at zero frequency is 
estimated by a linear extrapolation toward zero from the 
first two Matsubara frequencies. In Fig. 01 we present 
A{k,uj = 0) for two metallic cases U = 2.0,2.25. The 
Fermi surface is only weakly renormalized compared to 
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FIG. 3: Double occupancy as a function of U. The 
curves correspond (from bottom to top) to T/D — 
1/20,1/30,1/40,1/44. 



the U — value. Furthermore, the right part of Fig. ^ 
shows that as the Mott transition is approached the re- 
gions along the diagonal are characterized by a larger 
spectral intensity, thus a smaller electronic lifetime, than 
along the (0, tt) directions. Note that the existence of hot 
and cold spots along the Fermi surface have been the sub- 
ject of intensive discussions in the cuprate literature jl.lj. 
The mechanism for the momentum differentiation in our 
CDMFT study is as follows. When U increases, three sets 
of curves can be distinguished: for [/ = 2, 2.1, 2.2 (Green 
curves on Fig. I2I5|I . the on-site lattice self-energy ex- 
hibits metallic behavior, while the off-diagonal elements 
are very weak, leading to a DMFT-like result and to 
the A{k,uj = 0) plotted on the left in Fig. |3| For 
U = 2.25,2.29,2.31 (Blue curves). En is still metallic- 
like but the E"4,E23 are not negligible (see Fig. [SJ. We 
will refer to these points as the "intermediate region". 
For the largest U (2.34, 2.43, 2.5, red curves), the system 
is a finite temperature Mott insulator. In this regime 
the real part of E14 and E23 become large and have a 
sign such that it tends to restore the square symmetry 
of the lattice. In the intermediate region, E'/4 and S23 
become large and of comparable magnitude. This is a re- 
markable and unexpected effect, since the bare hopping 
amplitude between the 1,4 and the 2,3 sites are very dif- 
ferent (0. and .25 ) in our model. Hence, the anisotropy 
of the electron lifetime stems from cos(fca;) cos(fcj^). 
This quantity is maximum and negative around (0, tt) 
and (tt, 0) (hot spots) and is a minimum, i.e. zero, at the 
diagonal (7r/2,7r/2) (cold spots) and leads to the curves 
plotted on the right in Fig. ^ Note however that this 
is an intermediate temperature effect. Indeed a scenario 
where, at zero temperature and w — > 0, E'^ — > and 
E"4 ^ const (this would lead to a non Fermi liquid be- 
havior more pronounced around (0, tt) and (7r,0)) is for- 
bidden by causality (and thus by CDMFT): E" has to be 




FIG. 4: Distribution of low energy spectral weight in fc— space 
A{k,uj = 0) for T/D = 1/44, U/D = 2.0 (left panel) and 
U/D = 2.25 (right panel). The top panel are color plots to 
see the Fermi surface (red is high, blue is small) and bottom 
panel are 3d plots to see the variation of A. For intermediate 
U, cold and hot regions are visible around (f , f ) and (7r,0) 
respectively. 

negative as a matrix for all w, but if En(ix' = 0) = and 
E"4(a; = 0) < then the submatrix E^ with i,j e {1,4} 
would not be negative. 




FIG. 5: Imaginary part of E14 for T/D = 1/44 and U = 
2,2.1,2.2 (green curves), U = 2.25,2.29,2.31 (blue curves), 
U = 2.34, 2.43, 2.5 (red curves) and T/D = 1/44. 

More refined quantities are the variation of the quasi- 
particle residue and the effective mass on the Fermi sur- 
face. We find that quasi particle residue decreases with 
U, roughly linearly and becomes small near U k, Uc s& in 
DMFT. However, since it involves the derivatives of the 
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analytically continued self-energy, it can not be evaluated 
within our QMC results with enough precision to differ- 
entiate the behavior of Z in the hot and cold region. On 
the other hand, an interesting physical quantity, that can 
be reliably extract from the real part of the self energy 
at zero frequency, is: 

^ , dS'(fc,0) idEk 

where m* and m are respectively the renormalized and 
the bare effective mass, fc is along the Fermi surface and 
-^-^ means the derivative perpendicular to the Fermi sur- 
face, is plotted in Fig. Elfor fiD = 44 as a function 
of U for the hot and cold points. Note that — Z is al- 
ways equal to one in single site DMFT because of the fc- 
independence of the self energy. Instead within CDMFT 
we find that it decreases as we approach the Mott tran- 
sition mainly because of the increasing in modulus of the 
real part of E12 at zero frequency, see inset of Fig. 
This leads to a contribution to the /c-dcpcndent renor- 
malization of the effective mass which tends to counter- 
balance the increase of XjZ. Whether the effective mass 
at low temperatures remains finite as the Mott point 
is approached as in the large TV studies of Ref 14*, or 
whether it diverges in an albeit weaker fashion than in 
DMFT, remains an open question which will require a 
lower temperature study. 

Finally, analyzing the density of states using the Max- 
imum Entropy [l^ analytic continuation procedure we 
have found that the quasi particle peak found in DMFT 
in the metallic region is split near the Mott transition 
into two peaks, giving rise to a pseudo-gap |l6|| . 




FIG. 6: as a function of U at /3D = 44. Top (bottom) 

curve corresponds to the hot (cold) point. Inset :E'i2 vs. lOn 
for U = 2,2.1,2.2,2.25,2.29,2.31 (from top to bottom) and 
TjD = 1/44. 

To summarize, we found that the Mott transition, pre- 
vented by AF order in the unfrustrated model, is present 
in the frustrated case. The single site DMFT scenario 
characterized by a finite temperature Mott endpoint is 



compatible with our data. At low temperature a sub- 
stantial fc-dependence of the self-energy shows up and 
regions with significant variation of the electron lifetime 
as a function of k appear. Our analysis has been per- 
formed on a model relevant to the quasi-two dimensional 
materials of the k family [tR H^ . Indeed in Ref 
the Mott transition has been recently driven by pressure. 
It would be very interesting to carry out angle resolved 
photoemission experiments in these materials to verify 
our theoretical predictions, regarding the strong momen- 
tum dependence of the lifetimes. We conjecture, that 
this differentiation in momentum space is more general 
than the model we studied, and is direct consequence of 
the proximity to the Mott transition. Indeed it would 
be interesting to perform a similar study in a doped and 
isotropically frustrated Hubbard model relevant for high- 
Tc superconductors 0. 

The computations were carried out at the parallel clus- 
ter of the Physics Department at Rutgers. We thank 
C. Uebing for technical help with the computers, A. Per- 
ali for his participation in the early stages of this work. 
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